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Abstract 

The spin-network quantum simulator model, which essentially en- 
■ codes the (quantum deformed) SU (2) Racah-Wigner tensor algebra, 

is particularly suitable to address problems arising in low dimensional 
qh topology and group theory. In this combinatorial framework we im- 

plement families of finite-states and discrete-time quantum automata 
capable of accepting the language generated by the braid group, and 
whose transition amplitudes are colored Jones polynomials. The au- 
tomaton calculation of the polynomial of (the plat closure of) a link 
L on 2N strands at any fixed root of unity is shown to be bounded 
from above by a linear function of the number of crossings of the link, 
on the one hand, and polynomially bounded in terms of the braid 
index 2N, on the other. The growth rate of the time complexity func- 
tion in terms of the integer k appearing in the root of unity q can 
be estimated to be (polynomially) bounded by resorting to the field 
theoretical background given by the Chern-Simons theory. 

Key words: link invariants; braid group representations; Chern-Simons theory; 
quantum automata; Racah-Wigner algebra; spin-network simulator; topological 
quantum computation; U q (su(2)) representation theory. 
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1 Introduction 



The spin-network quantum simulator model |2] represents a bridge be- 
tween circuit schemes for standard quantum computation and approaches 
based on notions from Topological Quantum Field Theories (TQFT) |SHHE]- 
The spin-network computational space, naturally modelled as a graph for any 
fixed number of incoming spins, supports computing processes represented 
by families of paths and provides, on the one hand, a consistent discretized 
version of the topological quantum computation approach. On the other 
hand, such a quantum combinatorial scheme, which essentially encodes the 
(quantum deformed) SU(2) Racah-Wigner tensor algebra, turns out to be 
particularly suitable to address problems arising in (low dimensional) topol- 
ogy and group theory. The guiding idea of this paper is that the exponential 
efficiency that quantum algorithms may achieve with respect to classical ones 
proves to be especially relevant in problems in which the space of solutions is 
characterized by a structure definable in terms of the grammar and the syn- 
tax of a language, rather than algebraic or number-theoretic in nature. The 
spin-network setting provides a 'natural encoding' for classes of problems 
which basically share the combinatorial structure of the language underlying 
the (re)coupling theory of SU(2) angular momenta [§]. 

On the other hand, the Jones polynomial [7j is no doubt the most famous 
knot invariant in topology, a knot invariant being a function on knots (or 
links, namely circles embedded in 3-space) which is invariant under isotopy 
(smooth deformations) of the knot. Among its many connections to various 
mathematical and physical areas (see e.g. jH] for applications in statistical 
mechanics), we are mainly interested here in its relations with TQFT In 
the seminal paper PU| > Witten put link invariants in a field theoretical setting, 
showing that Jones polynomials arise as vacuum expectation values of Wilson 
loop operators in a three dimensional SU(2) Chern-Simons (topological) 
quantum field theory where the fundamental representation of the gauge 
group SU(2) lives on each component of the link. Such an invariant was 
extended to arbitrary representations living on the link components and in 
this paper we shall deal with such generalizations, referred to as 'extended' 
or 'colored' Jones polynomials [TT) I12j. 

From the (classical) computational side, it was proved that the exact 
evaluation of the Jones polynomial of a link L, V(L, uj) at uj = root of unity, 
can be performed in polynomial time in terms of the number of crossings of 
the planar diagram of L if u is a 2nd, 3rd, 4th, 6th root of unity. Otherwise, 
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the problem is #P-hard ^3] (the computational complexity class #P-hard 
is the enumerative analog of the NP class) . However, Kitaev, Larsen, Freed- 
man and Wang [3] showed that their 'topological' quantum computation 
setting, relying on the same TQFT quoted before, implicitly provides an 
efficient quantum algorithm for the approximation of the Jones polynomial 
at a fifth root of unity Unfortunately, this important algorithm was never 
explicitly formulated. This is particularly unfortunate since it is known that 
the approximation problem is BQP hard, and a quantum algorithm for this 
problem is thus of particular importance. 

Let us point out that recently Aharonov, Jones and Landau proposed 
an efficient quantum algorithm that approximates the problem of evaluating 
the Jones polynomial based, rather than on physical results from TQFT, on 
the path model representation of the braid group and the uniqueness of the 
Markov trace for the Temperley-Lieb algebra J3]- The argument is that 
the #P-hardness of the problem does not rule out the possibility of good 
approximations, and indeed these authors provide an efficient, explicit and 
simple quantum algorithm to approximate the Jones polynomial at all roots 
of unity for both the trace and the plat closures of a braid. 
Our strategy is quite different from theirs, since we shall basically provide a 
quantum (automaton) system whose internal evolution can be controlled in 
such a way that its probability amplitude gives the desired polynomial. 

As mentioned, one of the features of the Jones polynomial that will be 
used extensively is that it can also be defined via braids (a geometric N- 
braid is a set of iV strands with fixed endpoints in the plane). A braid can 
be 'closed up' to form a link by tying its ends together. In this paper we 
shall be interested in one of the two ways to perform such closures, namely 
the plat closure of the braid, and hence consider extended Jones polynomials 
associated with such link diagrams, cfr. Fig. ^ 

On a broader front, the study of braid groups and their applications is a 
field which has attracted great interest from physicists, mathematicians and 
computer scientists alike (cfr. jT3] for an updated review). Besides for its 
value in studying the braids in a theoretical framework, applications to knot 
theory have been known for years, while applicability to the field of cryptog- 
raphy has been realized recently The analysis of algorithmic problems 
related to braid group has thus acquired a great practical significance, in 
addition to its intrinsic theoretical interest. 

The approach we present here exploits a g-braided version of the original 
spin-network setting [2] to make it accept the language of the braid group 
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Figure 1: A plat presentation of the borromean link. 



and to deal with link polynomials (see for a presentation of some pre- 
liminary results). As pointed out before, the 'physical' background provided 
by the 3D quantum 577(2) Chern-Simons field theory plays a prominent 
role, because our computational scheme is actually designed as a discretized 
conterpart of the topological quantum computation setting proposed in [H]. 
Moreover, this framework is exactly what is needed to deal with (normalized) 
577 (2)-colored link polynomials expressed as vacuuum expectation values of 
composite Wilson loop operators, on the one hand, and with unitary rep- 
resentations of the braid group, on the other. These expectation values, 
in turn, will provide a bridge between the theory of formal languages and 
quantum computation, once more having as natural arena for discussion the 
g-braided spin-network environment. We are going to implement families of 
finite states (and discrete time)-quantum automata capable of accepting the 
language generated by the braid group, and whose transition amplitudes are 
colored Jones polynomials. More precisely, our results will be interpreted in 
terms of 'processing of words' -written in the alphabet given by the gener- 
ators of the braid group- on a quantum automaton in such a way that the 
expectation value associated with the internal automaton 'evolution' is ex- 
actly the extended Jones polynomial. The quantum automaton in question 
will in turn correspond to a path in the g-braided spin-network computa- 
tional graph. The calculation of the polynomial of (the plat closure of) a 
link L on 2N strands will be shown to be bounded from above by a linear 
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function of the number of crossings of the link, on the one hand, and poly- 
nomially bounded in terms of the braid index 2N ', on the other. Notice that 
the growth rate of the time complexity function in terms of the integer k 
appearing in the root of unity q can be easily estimated to be (polynomially) 
bounded by resorting to the TQFT background, since k is nothing but the 
Chern-Simons coupling constant. 

We shall leave as open problems the analysis of the complexity of the 
preparation of (initial and final) states as well as the efficient implementation 
of the individual automaton transition functions, which might be addressed 
by means of approximating (classical or quantum) algorithms. 
In conclusion, we argue that our field theoretical approach could be further 
generalized, by suitable modifications of the braiding prescriptions in the 
spin-network scheme, to deal with 2-variables link polynomials such as the 
HOMFLY invariant |18j . related to the partition function of Potts model jH]. 

The content of the paper is, as far as possible, self contained. In section 2 
we briefly recall the definitions of classical and quantum languages and finite 
states-automata. In section 3 we give a review of the spin-network computa- 
tional framework modelled on the Racah-Wigner tensor algebra of SU (2). In 
section 4 we deal with the g-braided version of the spin-network simulator, 
which relies on the tensor algebra of U q (su(2)) (at q = root of unity). Section 
5 is splitted into two parts: in 5.1 we review the 'quantum group approach' 
(and related i?-matrix) to the study of (unitary) braid group representations 
and 'quantum' link invariants; in 5.2 we present the field-theoretical back- 
ground (Chern-Simons TQFT, Wess-Zumino boundary theory, composite 
Wilson loop operators and their expectation values) trying to resort to geo- 
metric intuition rather than to a deep knowledge of techniques in quantum 
field theory. In section 6 we explain in details the automaton calculation of 
the extended Jones polynomial. 

2 Classical and quantum formal languages 

The theory of automata and formal languages addresses in a rigorous way 
the notions of computing machines and computational processes. If A is an 
alphabet, made of letters, digits or other symbols, and A* denotes the set of 
all finite sequences of words over A, a language C over A is a subset of A*. 
The length of the word w is denoted by \w\ and Wi is its i'th symbol. The 
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empty word is and the concatenation of two words u, v is denoted sim- 
ply by uv. In the sixties Noam Chomsky introduced a four level-hierarchy 
describing formal languages according to their structure (grammar and syn- 
tax): regular languages, context-free languages, context-sensitive languages 
and recursively enumerable languages. The processing of each language is in- 
herently related to a particular computing model (see e.g. [HI] fo r an account 
on formal languages). Here we are interested in finite states-automata, the 
machines able to accept regular languages. 

A deterministic finite state automaton (DFA) consists of a finite set of 
states S, an input alphabet A, a transition function F : S x A — > S, an 
initial state s in and a set of accepted states S acc C S. The automaton starts 
in Si n and reads an input word w from left to right. At the z— th step, if 
the automaton reads the word Wi, then it updates its state to s' = F(s,Wi), 
where s is the state of the automaton reading w%. One says that the word 
has been accepted if the final state reached after reading w is in S acc . 
In the case of a non-deterministic finite state automaton (NFA), the tran- 
sition function is defined as a map F : S x A — > P(S), where P(S) is the 
power set of S. After reading a particular symbol, the transition can lead 
to different states, according to some assigned probability distribution . If a 
NFA has n states, for each symbol a G A there is an n x n transition matrix 
M a for which (M a )ij = 1 if and only if the transition from the state i to the 
state j is allowed once the symbol a has been read. 

Generally speaking, quantum finite states-automata (QFA) are obtained 
from their classical probabilistic counterparts by moving from the notion 
of (classical) probability associated with transitions to quantum probability 
amplitudes. Computation takes place inside a suitable Hilbert space through 
unitary matrices and a number of different models have been proposed, see 
e.g. [2niE]; just to mention a couple of them. Following [21], the measure- 
once quantum automaton is a 5-tuple M = (Q, E, S, qo, q/), where Q is a 
finite set of states, E is a finite input alphabet with an end-marker symbol 
# and 5 : Q x E — > Q is the transition function. Here <5(q, a, q') is the 
probability amplitude for the transition from the state q to the state q' 
upon reading the symbol a. The state qo is the initial configuration of 
the system, and q/ is an accepted final states. For all states and symbols 
the function 5 must be unitary. The end-marker # is the last symbol of 
each input and computation terminates after reading it. At the end of the 
computation the automaton measures its configuration: if it is an accepted 
state then the input is accepted, otherwise is rejected. The configuration of 
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the automaton is in general a superposition of states in the Hilbert space 
where the automaton lives. The transition function is represented by a set 
of unitary matrices U a (a G £), where U a represents the unitary transition 
of the automaton reading the symbol o . The probability amplitude for the 
automaton of accepting the string w is given by 

}m (w) = (q/| U w |q ) , (1) 

and the explicit form of fuiw) defines the language C accepted by that 
particular automaton. If P denotes the projector over the accepted states, 
the probability for the automaton of accepting the string w is given by 

p M (w) = \\P\ij w )\\ 2 (2) 

where \ip w ) = U w |q ). 

3 The quantum spin— network simulator 

The spin network model of computation was introduced in P and worked 
out in [2] as a general framework for processing information in the quan- 
tum context and is essentially modelled on the combinatorics of the Racah- 
Wigner algebra of SU(2). The spin-network can be seen as a collection of 
graphs & n (V,E) parametrized by an integer n (n > 2), where n + 1 is the 
number of incoming angular momentum variables, each associated with an 
irreducible representation (irrep) of SU(2), {ji} G {0,1/2,1,3/2...} in h 
units (we choose units in which H — 1). On the physical side, these n + 1 
basic variables enter in the construction of different sets of (pure angular 
momenta) eigenspaces selected according to the different types of quantum 
interactions we whish to simulate. The fact that physical interactions in many 
(conservative) quantum systems can be well modelled on (combinations of) 
two-body interactions 22. opens the possibility of calling into play the pow- 
erful algebraic-combinatorial setting underlying SU(2) binary coupling and 
recoupling theory (cfr. [B] and the original references therein). 
Before going into some more details on this realization of the spin-network 
graphs, let us point out that the combinatorial structure encoded into the 
Racah- Wigner algebra is actually shared by other discrete structures. 

A first type of realization is purely graph-theoretical. The vertex set V 
of the graph (5 n (V, E) can be identified with the set of (rooted) binary trees 
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with ri + 1 labelled leaves where the leaves (terminal nodes) and the internal 
nodes of the trees are labelled by integers and half-integers G |N, cfr. Fig. 
121 Undirected edges between vertices are drawn whenever a pair of vertices 
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Figure 2: Two labelled binary trees on (n + 1) = 4 leaves. Such trees are in 
one-to-one correspondence with the vertex set V of the graph <8>3(V, E) 

(labelled trees) are connected by two kinds of topological elementary moves, 
namely twist and rotation, illustrated in Fig. El 

The resulting graph, known as Twist-Rotation graph, is depicted for 
n + 1 = 4 in Fig. 0] and its combinatorial properties are analyzed in [23] and 
in Appendix A of |2j. 

Another realization of the spin-network is in terms of words endowed with 
pairs of parentheses representing a non-commutative and non-associative bi- 
nary operation. In this case the vertices of the graph <3 n (V, E) are associated 
with words w made of letters from the alphabet { |N U pairs of labelled 
parentheses (-) }, e.g. 

w = {{{h,32)ki,h) k2 \3uh G | N with jx + j 2 + . .. j n+1 = J, (3) 

where J is the label assigned to the root. Two vertices are connected by an 
edge if it is possible to switch from one to the other either by swapping the 
elements inside a parenthesis, (a, b) c «-> (b, a) c , or by changing the parenthe- 
sization structure ((■ , -) fel , -) fc2 <-►(•,(•, ■)h 1 ) h2 - 

Coming back to the Racah-Wigner setting, the interpretation of the spin- 
network graph goes on as follows. There exists a one-to-one correspondence 
{v (b)} < — ► {H J n (b)} between the vertices of E) and the computational 

Hilbert spaces of the simulator. The label b has the following meaning: 
for any given pair (n, J), all binary coupling schemes of the n + 1 angular 
momenta {J^}, identified by the quantum numbers ji, ■ ■ ■ , j n +i (summing 
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rotation 



Figure 3: A twist corresponds to the interchange of either two leaves or two 
subtrees (top). A rotation consists in a change of the coupling scheme of 
either three leaves or subtrees (bottom). 

up to a total J) plus ki, . . . , fc n _i (corresponding to the n — 1 intermediate 
angular momenta {Kj}) and by the brackets defining the binary couplings, 
provide the 'alphabet' in which quantum information is encoded (the rules 
and constraints of bracketing are instead part of the 'syntax' of the resulting 
coding language). The Hilbert spaces TL^ (b) thus generated are spanned by 
complete orthonormal sets of states with suitable quantum number label set 
such as, e.g. for n = 3, { ((h(hh) kl ) J±) J ' ((^2) ^(jsk) k ,) j}- 

More precisely, for a given value of n, H^(b) is the simultaneous eigenspace 
of the squares of 2(n+l) Hermitean, mutually commuting angular momentum 
operators Ji, J 2 , J3, ...,J n+ i with fixed sum Ji + J 2 + J3 +. . .+J n+ i = J, 
of the intermediate angular momentum operators K 1; K 2 , K 3 , . . . , K n _! and 
of the operator J z (the projection of the total angular momentum J along 
the quantization axis). The associated quantum numbers are ji , j 2 , • • • , jn+i] 
J; ki, /c 2 , . . • , k n _i and M, where — J < M < +J in integer steps. 
If W 1 ® W 2 <g> ■ ■ ■ ®W n ®W n+1 denotes the factorized Hilbert space, namely 
the (n+l)-fold tensor product of the individual eigenspaces of the (J^) 2 's, the 
operators Kj's represent intermediate angular momenta generated, through 
Clebsch-Gordan series, whenever a pair of J/s are coupled. As an exam- 
ple, by coupling sequentially the J/s according to the scheme (• • • ((Ji + 
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Figure 4: A portion of the Twist-Rotation graph (5 3 (V, i?) where only 30 
out of 60 vertices are shown (the picture can be completed by taking the 
mirror image of each tree at the antipodal vertex). The remaining 60 vertices 
are arranged into an isomorphic graph obtained by swapping one pair of 
labels, e.g. (a, b) — > (6, a). Solid edges represent rotations and dashed edges 
represent twists. 

J2) + J3) H 1- Jn+i) = J which generates (Ji + J 2 ) = Ki, (K x + J 3 ) = 

K 2 , and so on - we should get a binary bracketing structure of the type 
(• • • (((W 1 ® H j2 ) kl ®W' A ) k2 ® • • • <g> n jn+1 ) kn _ 1 ) j, where for completeness we 
add an overall bracket labelled by the quantum number of the total angu- 
lar momentum J. Note that, as far as j/s quantum numbers are involved, 
any value belonging to {0,1/2,1,3/2,...} is allowed, while the ranges of 
the fcj's are suitably constrained by Clebsch-Gordan decompositions (e.g. if 
(Ji + J 2 ) = Ki =>• |ji - j 2 \ < h < ji + 32). 

We denote a binary coupled basis of (n + 1) angular momenta in the JM- 
representation and the corresponding Hilbert space as 

{ I [Ji, 32, h, ■ ■ ■ , J„+i] b ; k\ , k\ , . . . , k h n _ x ;JM),-J<M<J} 

= H J n (b) = span{|b;JM) n }, (4) 
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where the string inside [ji, j 2 , j 3 , . . . ,j n +i] b is not necessarily an ordered one, 
b indicates the current binary bracketing structure and the k^s are uniquely 
associated with the chain of pairwise couplings selected by b. 
For a given value of J each TC^(b) has dimension (2J+ 1) over C, but Hilbert 
spaces corresponding to different bracketing schemes, although isomorphic, 
are not identical. They actually correspond to (partially) different complete 
sets of physical observables, namely for instance {J 2 , J 2 ,, J 2 2 , J3? J 2 , J z } 
and {J 2 , J|, J3, J 2 3 , J 2 , J z } respectively (in particular, J 2 2 and J23 cannot 
be measured simultaneously). On the mathematical side this remark reflects 
the fact that the tensor product ® is an associative operation only up to 
isomorphisms. 

For what concerns unitary operations acting on the computational Hilbert 
spaces (j3J), we shall consider here unitary transformations associated with 
recoupling coefficients (3nj symbols) of SU(2), thought of as j-gates in the 
present quantum computing context. As shown in [Bj, any such coefficient can 
be splitted into 'elementary' j-gates, namely Racah and phase transforms. 
A Racah transform applied to a basis vector is defined formally as 

U :\...{{ab) d c) f ...-JM) ^\...(a(bc) e ) f ...;JM), (5) 

where Latin letters a,b,c, . . . are used here to denote generic, both incoming 
(je 's in the previous notation) and intermediate (fcj 's) spin quantum numbers 
(this operation corresponds to a rotation in the Twist-Rotation graph, erf. 
Fig. El bottom and Fig. HJ). Its explicit expression reads 

\(a(bc) e ) f ;M) 



= ^ ( _ ir+b+c+/[(2rf+1)(2e + 1)] i/ 2 | « b d J | ((a6)dc)/;M)) (6) 

where there appears the 6j symbol of SU (2) and / plays the role of the total 
angular momentum quantum number. Note that, according to the Wigner- 
Eckart theorem, the quantum number M (as well as the angular part of wave 
functions) is not altered by such transformations, and that the same happens 
with any 3nj symbol. On the other hand, the effect of a phase transform 
$ (a twist operation on the Twist-Rotation graph, see Fig. El top and Fig. 
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EJ) amounts to introducing a suitable phase whenever two spin labels are 
swapped 

\...(ab) c ...;JM) = (-l) a+b - c \...(ba) c ...;JM). (7) 

These unitary operations are combinatorially encoded into the edge set E = 
{e} of the graph <3 n (V, E): E is just the subset of the Cartesian product 
(V x V) selected by the action of these unitary j-gates. 

In the framework described above, a computation is represented in a nat- 
ural way by a collection of step-by-step transition rules (gates), namely a 
family of 'elementary unitary operations' and we assume that it takes one unit 
of the intrinsic discrete time variable to perform anyone of them. Such pre- 
scriptions amount to select (families of) 'directed paths' in the spin-network 
computational space <5 n (V,E) x C 2J+1 , all starting from the same input 
state and ending in an admissible output state. A single path in the given 
family can be interpreted as a (finite-states) quantum automaton calcula- 
tion, once we select a particular encoding scheme for the problem we wish to 
address. 

By a directed path V with fixed endpoints we mean a (time) ordered 
sequence 

l v in )» = l v o )n -> |vi } n -> ► |v s ) n -> ► |v L ) n = |v out )„ , (8) 

where we use the shorthand notation |v s ) n for computational states (which 
are vectors expressed in the bases (J1J) and s — 0,1,2, ... , L(P) is the lexico- 
graphical labelling of the states along the path. Finally, L(P) is the length 
of the path V and L('P) • r = T is the time required to perform the process 
in terms of the discrete time unit r. 

A computation consists in evaluating the expectation value of the unitary 
operator U-p associated with the path V, namely 

( v out I U p I v in )«• 

(9) 

By taking advantage of the possibility of decomposing XJ V uniquely into an 
ordered sequence of elementary gates, JSJ) becomes 

L-l 

( v out I u p I v in )n = lY[(vs+i\U s ,s+i\v s } n \ P (10) 

s=0 
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with L = L(P) for short. The symbol |_ \v denotes the ordered product 
along the path V and each elementary operation is rewritten as U s>s+ i (s = 
0, 1, 2, . . . L('P)) to stress its 'one-step' character. Such expectation values are 
particular instances of the general expression (0) for the quantum amplitude 
of a finite-states automaton, once a suitable language has been encoded into 
the computational space of the spin-network simulator. 

4 (/—braided computational space 

As we shall see in the following section, the basic ingredient for addressing 
link invariants arising in the context of Chern-Simons field theory is the 
'tensor structure' naturally associated with the representation ring of the Lie 
algebra of a simple compact group which plays the role of the gauge group 
of the theory. In the case of SU(2) this structure is provided by (tensor 
products of) Hilbert spaces supporting irreducible representations together 
with unitary morphisms between them: these are exactly the objects col- 
lected into the Racah-Wigner algebra discussed in section 3. However, when 
dealing with (planar diagrams of) links we shall also have to specify the eigen- 
values of the braiding matrix to be associated with the crossings of the links 
and this extension can be achieved by 'braiding' the Racah-Wigner tensor 
category. In the present context, it is natural to take advantage of quan- 
tum group techniques in order to 'split' any phase transform ((7|) by assigning 
different weights -depending on a deformation parameter q to be defined 
below- to right and left handed twists. From the combinatorial viewpoint, 
this generalization corresponds to replace the spin network computational 
space <8 n (V,E) x C 2J+1 with its g-braided counterpart (see Fig. |3J) 



where the (classical) 6j symbol in any Racah trasform (jUJ) will become q- 
deformed. 

The tensor category we are going to introduce is associated with the quan- 
tum group U q (su{2)) (q = root of unity), namely the universal enveloping 
algebra of SU (2) endowed with additional structures which make it a quasi- 
triangular quasi-Hopf-*algebra (see e.g. 24J and other references therein). 
U q (su (2)) is an associative algebra generated by elements J + , J_ and J z 
which satisfy the commutation relations 



((<& n (V,E) x C 2J+1 ) X Z 2 ) ? , 



(11) 



[J Z ,J ± ] = ±J ± ;[J + ,J^] = [2J Z ] 



g ' 



(12) 
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Figure 5: A portion of the g-braided Twist Rotation graph (<5 3 (V, E) x Z 2 ) g : 
with respect to the unbraided situation, each twist has been splitted. 



where the g-integer [n] 9 is defined as [n] g = (q n l 2 — q~ n / 2 ) / (q 1 ^ 2 — q^ 1 ^ 2 )- 
U q (su(2)) is a deformation of the universal enveloping algebra of the Lie 
algebra su(2) since in the limit q — > 1 the above relations reduce to the com- 
mutation relations for the su(2) generators. 

A Hopf algebra-structure can be introduced by defining the coproduct ho- 
momorphism 

A : U q {su{2)) -> U q {su{2)) ® U q (su(2)), 
acting on J + , J_ and J z according to 

A ( J ± ) = J± ® qi + q'i ® J ± ; 

A(J 2 ) = J 2 ® 1 + 1 ® J 2 . 

The tensor algebra associated with U q (su(2)) can be worked out in practice 
as in the case of su(2), so that we have Hilbert spaces supporting irreducible 
representations, g-Clebsch-Gordan coefficients, g-Racah coefficients and so 
on. The crucial difference consists in the fact that the irreps label set acquires 
a cut-off, namely each label must be chosen in the set {0,1/2,1,3/2,..., 
k — 2}, where the integer k is related to the deformation parameter q by q 
= exp(—2in/k). Denoting by TL^ 1 and W q 2 the Hilbert spaces supporting two 
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irreps ji,j 2 , their (truncated) tensor product can be decomposed according 
to the Clebsch-Gordan series 

min{j 1 +j 2 ,k-j 1 -j 2 } 

n®K = k- (13) 

j=b'i— J2I 

As happens in the classical case, the two bases associated with the eigenspaces 
involved in the tensor product (fT3|) can be connected by means of Clebsch- 
Gordan coefficients according to 

\j m ) q = (^hmim 2 I jm) q \j x m 1 ) q \j 2 m 2 ) q , (14) 

mi,ni2 

where —ji < mi < ji, — j 2 < m 2 < j 2 , m = m\ + m 2 . Quantum CG 
coefficients ( ) q can be suitable normalized and satisfy orthogonality relations 

The quantum Racah transformation comes out when we consider different 
binary couplings in the tensor product Wg 1 <S> W q 2 ® W q a of three irreducible 
representations, as done in the classical case (c/r. © and ©). For instance 



= W 1 Ulirihl J12J23) ( [2jl2 + I]j2j23 + l]g ) V2 |j'l(j2j3)i 23 ; 3™) 

323 

(15) 

where there appear the g-dimensions of the irreps involved. The components 
of W q are the Racah coefficients of the algebra U q {su 2 ) and these symbols 
satisfy orthogonality relations, symmetry properties and identities which look 
like suitable g-deformations of the corresponding classical ones (and reduce 
to them in the limit q — > 1) The quantum Racah coefficient and the 

g-counterpart of the Wigner 6j symbol differ as usual by a phase factor, 
namely 



h 32 J12I 

33 3 323 J 



q 



Finally, we introduce the (differently normalized) symbol 

3jl 2 3j23 



i 3 ' 


32 


312 \ 


V 33 


3 


323 J 



[2J12 + 1] [2j 23 + 11 



(16) 
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which, on the one hand, enhances the matrix character of the quantum Racah 
transform (|15|) and, on the other, is particularly suitable to be generalized 
to deal with more than three incoming spin labels. 

5 Quantum invariants of links and unitary 
representations of the braid group 

Link invariants are functions on links (collections of knots, namely closed 
circles in 3-space) which depend only on the isotopy class of the link. An 
(ambient) isotopy can be thought of as a continuous transformation per- 
formed on the link embedded in IR 3 which deformes at will the shape of the 
link without cuttings. Let us point out preliminarly that the link invariants of 
polynomial-type we are going to address here are 'universal' in the sense that 
historically distinct approaches (R-matrix representations obtained with the 
quantum group method, monodromy representations of the braid group in 
2D conformal field theories, the quasi tensor category approach by Drinfeld 
and the 3D quantum Chern-Simons theory, see e.g. |26| l2*7j for reviews) are 
indeed different aspects of the same underlying algebraic structure. We shall 
focus in particular on the Chern-Simons setting jTUj since, on the one hand, 
it embraces the universal structure of (unitary) braid group representations 
shared by all the models quoted above and, on the other, can be naturally en- 
coded into the (braided) spin-network computational scheme. The (colored) 
link polynomials arising from SU (2) quantum CS theory can be referred to 
as 'extended' Jones polynomials, since the Jones polynomial [7] is recovered 
by selecting the fundamental (j = |) representation of 577(2) on each of the 
link components (or on each strand of the associated braid). Moreover, the 
topological quantum field approach is inherently related to low dimensional 
geometry since, for instance, suitable combinations of these invariants can 
be interpreted as topological invariants of hyperbolic 3-manifolds, obtained 
by surgery along framed links in the 3-sphere fTJ ^] . 

Let us point out that the definitions of link polynomials from Hecke (or 
Temperley-Lieb) algebra realizations of the braid group -exploited in 
in the quantum computational context- can be derived quite easily in the 
framework we are adopting here, since it can be shown that the associated 
invariants do satisfy the linear skein relations which characterize such real- 
izations pu] . 



16 



Before addressing a full fledged approach to 3D Chern-Simons theory, we 
pause a little bit digressing on i?-matrix representations of the braid group 
arising from quantum groups. The associated invariants of knots and links 
are commonly refereed to as 'quantum' invariants, since they are quantities 
depending on the deformation parameter q of the 'quantum group' under 
consideration. 

5.1 The quantum group approach 

Let be a (semi) simple Lie algebra, U q (g) its universal enveloping algebra 
and V a finite dimensional (complex) vector space in the associated tensor 
algebra (the prototype is of course the unitary tensor algebra of U q (su(2)) 
described in details in section 4). The representation theory of any such 
quantum group is naturally endowed with an invertible linear operator, the 
so-called .R-matrix 

R:V ®V ® V, (17) 
which satisfies the quantum Yang-Baxter equation 

(R ® /)(/ <g> R)(R ®I) = (I®R)(R ® ® R), (18) 

where both sides of the above expression are to be understood as linear 
transformations V ®V ®V — > V ®V ®V . 

The general procedure for constructing quantum invariants of oriented 
knots (or links) presented as closure (or platting) of braids can be outlined 
as follows. Consider an oriented knot diagram (namely the projection of a 
knot with orientation onto a fixed plane) and insert an horizontal line as 
depicted in Fig. El 

To each intersection point between the line and the diagram we assign 
either the representation space V or its dual V*, depending on whether the 
portion of the knot nearby the intersection is oriented upwards or downwards. 
The whole configuration of such points on the line turns out to be associated 
with the tensor product of the individual vector spaces (orderered from left 
to right). The connection with braid groups comes out when we consider 
two parallel horizontal lines intersecting the knot diagram. More precisely, 
the portion of the knot diagram between a pair of horizontal lines represents 
the geometric realization of a braid b, which in turn is an element of the 
Artin braid group B„, for some suitable n. B n has n generators, denoted by 
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V 



V* V 



V* 



Figure 6: The oriented trefoil knot cut by an horizontal line. We associate 
with the ordered set of the intersection points (from left to right) the tensor 
product V <8> V* <S> V <S> V*, where each factor is chosen in order to comply 
with the diagram orientation. 

{o"i, cr 2 , ■ ■ . , cr n _i} plus the identity e, which satisfy the relations 
o, o\, (\, oi if |i-j|>l = 1,2, ...,n- 1) 

Oi a i+ i <Ji = a i+1 OiO i+ i (i = 1,2, ...,n-2). (19) 

An element of the braid group is a word in the standard generators of B n , 
e.g. b = (T^ 1 a 2 o-^ 1 a 2 crf o 2 x <j\o 2 2 G B 4 ; the length |6| of the word b is 
the number of its letters. The group acts naturally on topological sets of n 
disjoint strands - ordered from left to right - in the sense that each generator 
CTj corresponds to the over-crossing of the ith strand on the (z + l)-th, and 
cr^ 1 represents the inverse operation (under-crossing) according to o~ x 
= o-rVj = e. 

On the other hand, when we represent B„ in the tensor algebra of U q (g), 
the action of a braid b is naturally associated with a linear operator T(b) 
connecting the vector spaces introduced above, see Fig. [7| 
Since T is a linear representation, we can simply specify its action on the 
standard generators {ai, cr 2 , . . . , cr n _i} to get {T(a 1 ),T(a 2 ), . . . , T(cr n _i)}, 
and extend this action to T(b) by linearity. The i?-matrix, namely the linear 
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Figure 7: The action of the braid group element b is represented as a map 
T{b) between the vector spaces living on the bottom and top lines. 

operator introduced in JT7J), is to be intended as the set of (elementary) cross- 
ing operators in some given representation T, constrained by the quantum 
Yang-Baxter relation (|18|). 

Knot theory is closely related to (representations of) braid groups owing 
to Alexander's theorem |28| . which states that every knot (or link ) L in the 
3-sphere S 3 = R 3 U oo can be presented (not uniquely) as a closed braid for 
some suitable n (to get a knot from the open braid of Fig. [7| we have to 
connect with arcs the lower and upper endpoints of each strand). We might 
also consider the plat presentation of a knot (characterized by the fact that 
the braid involved must possess an even number of strands), which is exactly 
the type of presentation depicted in Fig. El for the trefoil knot (see also Fig. 
^and Fig. ^2 hi the Appendix). Anyway, we can generate invariants of knots 
(links) for both types of presentations by taking some 'trace' of the operator 
T(b), where b is the braid associated with the given knot or link. The fact 
that the resulting quantities must depend only on the isotopy type of the 
knot can be suitably translated into the braid group-setting by resorting to 
the notion of invariance under Markov moves, and thus we should actually 
speak of 'Markov traces' (cfr. ^5J|2n] for reviews on knot theory and braid 
group). 

Summing up, the quantum group approach provides a purely algebraic 
construction of link invariants as (Markov) traces of representation matrices 
of the braid group in the tensor algebra of U q (g). Such invariants are poly- 
nomials in the deformation parameter q and its inverse 1/q. In the case of 
U q (su(2)) (q a root of unity), the associated g-braided Racah-Wigner algebra 
(discussed in section 4) is naturally endowed with Hilbert spaces and uni- 
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tary operators, namely the ideal arena to address (quantum) computational 
problems concerning both link polynomials and braid group. 

5.2 The C hern— Simons field theory approach 

A topological quantum field theory (TQFT) is a particular type of gauge 
theory, namely a theory quantized through the (Euclidean) path integral 
prescription starting from a classical Yang-Mills action defined on a suitable 
.D-dimensional space(time). TQFT are characterized by observables (corre- 
lation functions) which depend only on the global features of the space on 
which these theories live, namely they are independent of any metric which 
may be used to define the underlying classical theory. The geometrical gen- 
erating functionals and correlation functions of such theories are computable 
by standard techniques in quantum field theory and provide novel represen- 
tations of certain global invariants (for D-manifolds and/or for particular 
submanifols embedded in the ambient space) which are of prime interest. 
Let us recall in brief the basic axioms for a unitary TQFT in D = 3 before 
going through the case which is of interest here, namely SU(2) Chern-Simons 
theory 0. 

Denote by Si and S 2 a pair of 2-dimensional manifolds and by M 3 a 
3-dimensional manifold with boundary DM 3 = Si U S 2 (all manifolds here 
are compact, smooth and oriented). A unitary 3-dimensional quantum field 
theory corresponds to the assignment of 

i) finite dimensional Hilbert spaces (endowed with non-degenerate bilinear 
forms) 7isi an d to Si and S 2 , respectively; 

ii) a map (technically, a functor) connecting such Hilbert spaces 

W a ^ W a (20) 

where A4 3 is a manifold which interpolates between Si (incoming boundary) 
and S 2 (outgoing boundary). Without entering into details concerning a few 
more axioms (diffeomorphism invariance, factorization etc.) we just recall 
that unitarity implies that 

iii) if S denotes the surface S with the opposite orientation, then 7ig = Ti,^, 
where * stands for complex conjugation; 

iv) the mappings (J2(J|) are unitary and Z[A4 3 ] = Z*LV1 3 ], where M 3 denote 
the manifold with the opposite orientation with respect to M 3 . 
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The classical SU (2) Chern-Simons action for the sphere S 3 (which is the 
simplest compact, oriented 3-manifold without boundary) is given by 

kS{A) = ^-J tr(AdA+^AAAAA) (21) 

where A is the connection 1-form with value in the Lie algebra su(2) of the 
gauge group, k is the coupling constant, d is the exterior differential, A is the 
wedge product of differential forms and the trace is taken over Lie algebra 
indices. The partition function of the quantum theory corresponds to the 
map (|2T)|) restricted to the case of empty boundaries and is obtained as a 
'path integral', namely by integrating the exponential of i times the classical 
action (}2*Tj) over the space of gauge-invariant flat SU(2) connections (the 
field variables) according to the formal expression 

Z cs [S 3 ; k] = J [DA] exp j ^ S cs (A) } (22) 

where the coupling constant k is constrained to be a positive integer by 
the gauge-invariant quantization procedure and is related to the deforma- 
tion parameter q (see below). The generating functional ()22|). written for a 
generic compact oriented 3-manifold A4 3 with dM 3 = 0, is a global invari- 
ant, namely depends only on the topological type [TO] . 

The extension of (|22j) to the case of a manifold with boundaries, DM 3 ^ 0, 
requires modifications of the classical action (}2"Tj) by suitable Wess-Zumino- 
type terms to be associated with each boundary component [3Uj ■ However, we 
do not need here the explicit expression of such boundary action since what 
we are interested in are expectation values of observables in the quantized 
field theory which will just require the knowledge of (vectors belonging to) the 
boundary Hilbert spaces, cfr. i) above. In particular, it turns out that the 
gauge-invariant observables in the quantum CS theory are expectation values 
of Wilson line operators associated with oriented knots (links) embedded in 
the 3-manifold (commonly referred to as Wilson loop operators). Knots and 
link are 'colored' with irreps of the gauge group SU(2), restricted to values 
ranging over {0, 1/2, 1, 3/2, . . . , k — 2}, where the integer k is related to the 
deformation parameter q by q = exp(— 2in/k) (see section 4 for details on 
the U q (su{2)) representation algebra). 

The Wilson loop operator associated with a knot K carrying a spin-j 
irreducible representation is defined as (the trace of) the holonomy of the 
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connection 1-form A evaluated along the closed loop K G S 3 , namely 



Wj [K] = t rj Pexp j> A, (23) 

where P is the path ordering. For a link L made of a collection of s knots 
{Ki}, each labelled by an irrep, the expression of the composite Wilson 
operator reads 

s 

W Jlj2 ... is [L] = J] W j; [iq. (24) 
1=1 

In the framework of the path integral quantization procedure, expectation 
values of observables are defined as functional averaging weighed with the 
exponential of the classical action. In particular, the functional average of 
the Wilson operator (J24*)l is 



n-Js L J frr. „ lfc 5cs (A) 



f[DA] 



e 4 ^ 



where Sqs (^4) is the CS action for the 3-sphere given in ()21|) and the gener- 
ating functional in the denominator will be normalized to 1 in what follows. 
It can be shown that this expectation value, which essentially 1 coincides with 
the extended (colored) Jones polynomial (H1H2], depends only on the iso- 
topy type of the oriented link L and on the set of irreps {ji, j s } (note also 
that S[L] — S[L], where L is obtained from L by reversing the orientation). 

The explicit evaluation of (|23j) can be carried out in several ways, by 
resorting to either field-theoretic methods, quantum group approaches (out- 
lined above) or through combinatorial state sum functionals. For future 
convenience we just sketch here the approach which relies on the extension 
of CS quantum theory -endowed with a Wess-Zumino conformal field theory 

1 These polynomials are actually invariants of 'regular' isotopy, which represents a re- 
stricted form of 'ambient' isotopy defined at the beginning of this section. The connection 
between £j 1 ...j s [L] and the genuine colored Jones polynomial is given by Jj 1 ...j B (L, q 1 ^ 2 ) = 
|q-3iu(L)/4y^i/2 _ g-i/2-jj £j i once suitable normalizations for the unknots have 

been chosen. Here w(L) is the writhe associated with the planar diagram D(L) of the link 
L, defined as w(L) = ^2 p s(p)- The summation runs over the self crossing points of D(L) 
and e(p) = ±1 according to simple combinatorial rules (see e.g. |23|)- The writhe is easily 
recognized from the link diagram by simple counting arguments, so that computational 
problems involving both link invariants belong to the same complexity class. 
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on its boundary- to the case in which the boundary components are inter- 
sected by knots or links, namely become 2-manifolds with punctures (note 
that this setting is closely related to the topological quantum computation 
approach j^j). The basic geometric ingredients can be easily visualized as 
in Fig. |H1 where a portion of a 3-dimensional manifold M. z (technically, a 
handlebody decomposition) is shown, together with an incoming boundary 
Si and an outgoing boundary S 2 made of two disjoint components, S 2 and 
S 2 . A portion of some knot (link) embedded in the ambient 3-manifold is 
also depicted, and its intersections with the boundaries are 'punctures' which 
inherit the irreps labels from the associated (Wilson) lines. 




Figure 8: A portion of an oriented 3-manifold with one incoming boundary 
and two outgoing boundaries. Lines belong to some knot (or link) embedded 
in the manifold and intersect the 2D boundaries in some points (punctures). 

According to the axioms of TQFT, we may associate with each boundary 
a (finite-dimensional) Hilbert space, that is Ti^ for the incoming boundary 
and H% 2 = ® ( nere > for simplicity, we do not explicitate the labels 
of puncures). The Chern-Simons unitary functional (see axiom ii)) is a state 
in the tensor product of these Hilbert spaces or, more precisely, 

Z cs [M 3 ;k} : W El ^Hx 2 

=> Z cs [M 3 ;k] e W El ®H*x 2 , (26) 

where in the last row we have used also axiom iii) since the incoming and 
outgoing boundaries must be endowed with opposite orientations. Moreover, 
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such type of expression is compatible with the quantum group approach 
outlined in section 5.1 since the the Chern-Simons mapping in (J26|) . when 
restricted to punctures, induces automatically (unitary) representations of 
the braid group in the tensor algebra of U g (su(2)). 

Finally, it can be shown [5| that the conformal blocks of the SU(2)e Wess- 
Zumino field theory living on the boundaries with punctures actually pro- 
vide the basis vectors for the Hilbert spaces introduced above (the level £ 
of the WZ model is related to the deformation parameter q according to 
q = exp{— 2m/ {£ + 2)}, and in turn £ is related to the coupling constant 
k(> 3) of the CS theory in the bulk by £ = k — 2). In the following section 
we shall carry on the explicit construction of such bases, which will allow us 
to recast the expectation value of the composite Wilson operator (|25j) into a 
form suitable to be handled for computational purposes. 

6 Automaton calculation of extended 
Jones polynomials 

As anticipated at the end of the previous section, we start with the con- 
struction of the basis vectors which will enter into the explicit expression of 
the expectation value of the composite Wilson operator (J25j) emerging from 
quantized 3D Chern-Simons theory. We use here the general setting given 
in [31] since it can be easily adapted to the g-braided spin-network scheme 
of section 4. 

Consider an oriented link L embedded in the 3-sphere, S 3 = 1R 3 U oo, en- 
dowed with a plat representation, namely presented as the closure of an 
oriented braid with 2N strands (cfr. Figg. HI lol and I11J1. If we remove two 
open three-balls from S 3 we get two boundaries, Si and S 2 , both topo- 
logically equivalent to S 2 , but with opposite orientations, (S ,3 ;Si,S 2 ) = 
(S 3 ; S 2 , S 2 ) (recall from section 5.2 that an SU(2)i Wess-Zumino conformal 
field theory is naturally associated with the oriented boundary surfaces). We 
can accomodate in such an ambient manifold, 2N 'unbraided' Wilson lines 
carrying irreps ji, ■ ■ ■ , 32N1 starting from the incoming (lower) boundary 
and ending into the outgoing (upper) one (punctures inherit the labellings 
ji G {0,l/2,l,3/2,...£} from the strands of the braid). Denote this 'identity' 
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colored oriented braid as 



Vi ( I 1 I 2 ' ' ' l 2N ) , (27) 

V h h ■■■ 32N J 

where ji = (ji,ei) i = 1,2, ...2N represents the spin ji together with an 
orientation = ±1 for a strand going into or away from the boundary, while 
stars over the symbols represent here the opposite choice of the orientation. 
In order to generate an arbitrary (oriented) braid i/g out of the identity 
braid uj we have to apply a braiding operator, denoted by the symbol B and 
written in terms of generators B\, B 2 , B 2 n-i to be defined below, starting 
from the lower boundary. With such prescription we shall get the braid 

"B {} f f f ), (28) 

where the labels have been ordered according to the requirement of having a 
plat presentation for the associated oriented link. 

Our goal will consist in recasting the expectation value of the composite 
Wilson operator, written in functional terms in (J25|) . into an expression which 
contains a quantity of the type 

(<P\B (f f f f (29) 

where B( ::: ) is the operator associated with the oriented braid (|28j). The 
shorthand notations \(f) > and \(f) > represent correlators of 2N primary fields 
in the SU (2)^ WZ theory, to be interpreted here as states belonging to the 
boundary Hilbert spaces associated, respectively, with the incoming and out- 
going Hilbert spaces of the underlying CS theory (cfr. the axioms for TQFT 
in section 5.2). 

The basis vectors to be associated with the incoming boundary can be de- 
noted in general as 

\\ji,32, J2jv-i, J2iv]; [k; h]; 0,0), (30) 

where the last two entries are the quantum numbers JM = (00) for a sin- 
glet state of the total angular momentum and the first string represents 
the incoming spin variables (we drop the hat on oriented objects whenever 
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not necessary). The second group of entries denotes particular sets of in- 
termediate angular momentum labels, arising from binary couplings of the 
2N primary fields j's, with a bipartite structure represented by the symbols 
k = fci, k%, . . . and h = hi, h 2 , . . ., chosen in order to comply with the rules 
described below. Before addressing the latter in general, let us have a look at 
the simplest non trivial case of N = 4 incoming spin labels. The underlying 
admissible combinatorial structures are depicted in Fig. El 




Figure 9: Combinatorial realization of the two basis sets in the case N = 4 
as labelled binary trees. They are connected by a duality matrix. 

In this case we have just two types of basis vectors, related by the so- 
called duality matrix of WZ theory 



Ibi,j2,j3,j4]; [M;-];00) = ^ 



1 ' U,j2,j 3 ,j4]; [&';§]; oo) 



k' 



Jl 32 

33 34 , „ 

(31) 

where the vectors on the left hand side do not contain any h-label, namely 
[k, h] = [k,k] — ], while on the right we have the combination [k', h'] = [k'; |]. 
The array ( :: | :) q in (|31|) is the (normalized) g-Racah symbol of U q (su(2)) 
introduced in (pH)j) of section 4 (with respect the notation used there, from 
now on we drop the subscript q on vectors). 

In the general case of 2N incoming spin labels, the two combinatorially 
distinct bases which have to be involved are specializations of the vectors 
in (J3U|) to the two configurations depicted in Fig. The extension of the 
duality transformation (J31j) to the case of an arbitrary (even) number of 
incoming spins can be done by resorting to two types of more complicated 
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Figure 10: Coupling binary trees representing the combinatorics of the two 
sets of bases in the case of a generic N, see (J35j) and 

arrays, which can be represented as 



3 1 

33 



32 

U 



J2N-5 32N-4 
32N-Z 32N-2 
\ 32N-X 32N 



h 
k 2 



h \ ( 

h 2 



J 



32 

H 



33 

3s 



32N~i 32N-3 
32N-2 32N-1 
\ 32N jl 



r 2 



8 1 

S-2 



\ 



r N - 2 S N -2 
r N -i stv_i 



/ 



(32) 

where the matrix indices -to be involved in summations whenever transfor- 
mations which generalize (jHljl are implemented- are listed in the right hand 
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side of the arrays. 

As happens in the standard Racah-Wigner setting, it can be shown that each 
of these arrays, which represent the g-deformed counterparts of SU(2) 3nj 
coefficients, can be decomposed in terms of g-Racah transformations (|31|) . 
cfr. sections 3,4 of and [3*T] . 

When the outgoing Hilbert space is considered (corresponding to the 
boundary S 2 endowed with the opposite orientation with respect to Ej) we 
have to introduce bra-type bases which are dual (and orthonormal) with 
respect to the bases in (J3*U|) . With an obvious choice of notations we set 

([ji, 32n\\ [k;h];0,0| [ji, j 2 jv]; [k'; h'] ; 0, 0) = 6 k> y5 hjW , (33) 

where, as before, h, k, h', k' represent multi-indices to be associated with the 
admissible configurations of binary coupled spins and there appear multiple 
Kronecker deltas. 

The discussion above was aimed to recognizing the crucial fact that the 
basis vectors 

I bi >h,-, J2N-1, 32n] ][k 1) ...,k N \h X: ... : h N _ 2 ] ; 0, 0) (34) 
are eigenfunctions of the odd braiding operators B 2 i-i, while the basis vectors 

\\ji,h, ...Jzn-i, jW]; [n, ...,rjv-i; s u s N -i]; 0, 0) (35) 

are eigenfunctions of the even braiding operators B 2 i- 

The explicit expressions of the eigenvalues, in the odd and even case respec- 
tively, read 

KChi-uhi) = Ai +) (J,J") = (-y+i'-* q^ +c i>y 2+c ^>)-°*/* for ee' = +1 

KCkuhi+i) = (X^UJ'T 1 = (-)^^'l-^lfe^'l/ 2 -^ 2 for ee' = -1. 

(36) 

Here / = 1, 2, . . . , N — 1, q is the deformation parameter, z G {ki, k 2 , ■ ■ ■ , kjy, 
r l5 r 2 , . . . , rjy-i}, c z = z{z + 1) is the quadratic Casimir for the spin-2; repre- 
sentation and e, e' denote the orientation of the strands labelled by j and j', 
respectively. Thus X^(j,f) is the eigenvalue of the matrix which performs 
a right handed half-twist in contiguous strands with the same orientation, 
while Xz (j, j') is the eigenvalue of the matrix which performs a right handed 
half-twist in strands with opposite orientation. 
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The explicit expression of the formal expectation value given in (|29|) above 
gives, after normalization according to the standard conventions (cfr. |31|). 
the extended Jones polynomial of the colored link L associated with the braid 
(J2Ej), namely 

N 

V jlj2 ... jN [L; q] J] [2ji + 1], x 
i=i 

([Ul,..J N J%];[0-,0]-A0\B^ d S'" ff)\[j l ,jl...J Nl j* N ]-,[0,0];0 1 0) 7 

'1 'i ■•• 'JV 'at 

(37) 

where [2ji + l] q is the g-dimension of the U q (su(2)) irrep j% defined in section 
4. The operator B{ ::: ) is expressed in terms of (a finite sequence of) the el- 
ementary braiding operators {B 2 i-u B21}, suitably changed into the current 
odd (even) basis by acting with a g-duality matrix, whenever an even (odd) 
vector of type (j33|) ((jMJ), respectively) is encountered. According to the ex- 
pressions (132)) and (J36)) for the admissible operations, the running variable 
of the polynomial is given by q = exp{— 27ri/(£ + 2)} for any integer £ > 3. 
Moreover, the above expectation value is to be interpreted as a trace over 
free spin labels. This feature derives of course from the geometric construc- 
tion of the plat presentation of the link L outlined at the beginning of this 
section, since the colored oriented braid (J2*S|) has to be 'closed up' to get the 
associated link. More precisely, corresponding strands must not only match 
pairwise with the correct orientations (as is made manifest by our notation), 
but j and /-type labels have to be appropriately identified (traced) in pairs, 
namely the Vs in (137)) are not new independent labels. 

Coming back to the 3-dimensional picture, such trace procedure amounts 
to gluing back the two opposite-oriented boundary 2-spheres, paying atten- 
tion to the coloring of the punctures, to end up with the same 3-sphere we 
started from. As already pointed out, the resulting link polynomial, arising 
as vacuum expectation value of the composite Wilson loop operator (j23|) in 
the quantum SU(2) CS theory for M. 3 = S 3 , is automatically an invariant 
of regular isotopy (cfr. the remarks in the footnote of section 5.2). 

This long technical discussion about the derivation of the extended Jones 
polynomial is nothing but the necessary premise to address the main issue of 
the present paper, namely the analysis of the connections among the theory 
of formal languages (section 2) and the spin-network computational scheme 
(sections 3 and 4), on the one hand, and braid group and links invariants on 
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the other. We are going to interpret the results established so far in terms of 
'processing of words', written in the alphabet given by the generators of the 
braid group, on a quantum automaton in such a way that the expectation 
value associated with the 'evolution' of the automaton is precisely the ex- 
tended Jones polynomial. The quantum automaton in question will in turn 
correspond to a path in the g-braided spin-network computational graph. 

In order to comply with the requirements for a finite-states quantum 
automaton described in section 2, we have to provide explicitly the 5-tuple 
(Q, E, S, qo, F), where F represents a set of acceptable final states. Now Q is 
a finite set of states belonging to the Hilbert spaces of the tensor algebra of 
U g (su(2)) described in section 5.1, whose combinatorial content was depicted 
in Fig. [TU| Labels of N irreps in this quantum group are associated with the 
strands of the plat presentation of the link L, and B2N is the braid group to 
be selected. £ is the alphabet made of the 2N — 1 generators of B 2 at: each 
generator (and its inverse) represents a letter of the alphabet, and words are 
written as composition of these elementary braids. The function 5 denotes 
a set of unitary matrices defining the transition rules and there is of course 
one matrix for each letter of the alphabet since we are linearly representing 
the braid group in the tensor algebra U q (su(2)). 

As initial state qo we pick up one particular binary-coupled state, namely 

|[0'i,i 2 ),-,0 2 7V-i,i 2 iv)];[0;0];0,0), (38) 

where, with respect to the generic expression for an odd basis vector given in 
pijl. we choose [k; h] = [0; 0], namely we select a 'multi-singlet' intermediate 
state. 

The set of final states F is constrained by the topological properties of the 
plat presentation, namely final states may differ from qo by a permutation 
on the string (jij2 • • ■ 32n)- Thus we can actually build up a family of N\ 
automata out of one initial state q . By acting with the symmetric group on 
the binary parenthesization structure of (}3*H|) we may get, for instance, the 
singlet final state 

|[((((0*3, h), UiJd), UwJb)), •••), (j2N-i,m ; [0, 0]; 0, 0) . (39) 

The unitary transition rules codified in the set S are: 

• if the automaton is in an even (odd) state and it reads an even (odd) 
braid generator, then the system evolves with the R-matrix associated 
to the proper braid generator, B 2 i-i or B 2 i, see ()34|) . (J35|) and ()36|): 
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• if the automaton is in an odd (even) state and it reads an even (odd) 
braid generator, then the system evolves with the proper duality trans- 
formation (see ([32))) to update the actual state into the configuration 
consistent with the parity of the given braid generator. As pointed out 
before, such a transformation can be splitted into a finite sequence of 
elementary duality (g-Racah) transformations of the type (|31|) . 

Once a final state q/ has been selected (the right permutation can be sin- 
gled out in a fast way even by a classical machine) the evaluation of the 
polynomial ()3*7jl is carried out by the automaton in a number of steps lin- 
ear in the length \B\ of the 'word' B. Since \B\ is the sum of the numbers 
of elementary braiding operators and g-duality transformations entering the 
explicit expression of B, the length of the word is bounded from above by 
a linear function of the number of crossings of the plat presentation of the 
associated link L. On the other hand, in the worst case we have to perform 
one duality transformation ()32|) before applying each elementary braiding 
operator Bi G {B 2 i-i, B 2 i}. This happens, for instance, in the evaluation of 
the Jones polynomial of the trefoil knot illustrated in the Appendix. This 
latter remark lead us to conclude that the time complexity function (more 
precisely, the number of computational steps in our automaton calculation) 
equals the length of B and is bounded according to 

\B\ < c(N)k(L), (40) 

where c(N) is a positive number depending on the number of generators of 
the braid group B 2 at and k(L) is the number of crossings in the plat presen- 
tation of the link L. 

We can estimate c(N) by observing that any such automaton can be uniquely 
associated with a particular path V in the g-braided spin-network compu- 
tational space (<b n {V,E) x C 2J+1 ) x Z 2 for n = 2N - 1. In particular, the 
maximum number of elementary g-Racah transforms entering a duality ma- 
trix of type ()32|) must coincide with the number of Racah transforms entering 
into one (classical) 3nj symbol since the combinatorics of such operations is 
manifestly the same. Hence we may exploit results from graph theory which 
tell us that the Rotation graph <5 n (V, E) -obtained from <8 n (V, E) by ignor- 
ing twists (or braidings)- has a diameter of the order n Inn (the diameter 
is defined as the maximum over the set of distances between pairs of ver- 
tices, where the distance is the minimum number of edges connecting two 
given vertices) (c/r. (221 and appendix A of [2 for a complete discussion 
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of the spin-network combinatorics). Clearly the 'distance' between the cur- 
rent basis and the eigenbasis with the right parity cannot exceed the above 
maximum distance, and consequently the factor c(N) in (J4*U|) grows as 

c(N) ~ (2JV - 1) In (2N - 1), (41) 

namely polynomially in the number of strands of the link. 

A deeper connection with the g-braided spin-network computational 
scheme comes out however when we recognize that the expectation value 
(J37j) representing the extended Jones polynomial is not only the quantum 
transition amplitude of a finite states-automaton, as pointed out before, 
but complies also with the expectation value (|lUj) to be associated with a 
path V in the g-version of the spin-network computational space. In this 
new perspective, what we are really doing is to 'encode' the combinatorial 
structure underlying quantum SU{2) Chern-Simons field theory (and the 
associated WZ boundary theory) at some fixed level I into the abstract q- 
braided S77(2)-spin-network for q = exp{— 2n/ (£ + 2)}. This does not mean, 
of course, that we have set up a quantum algorithm for the extended Jones 
polynomial in the strict sense, since the encoding map could not be 'effi- 
ciently' represented (nor efficiently approximated) with respect to standard 
models of computation (Boolean circuits, Turing machines). We provide, 
however, a quantum system whose evolution can be controlled in such a way 
that its probability amplitudes give the desired link polynomials. 

The crucial issue of constructing a bona fide quantum algorithm is under 
investigation. It will require in particular: i) (efficient) encoding schemes for 
binary coupled states; ii) (efficient) algorithms to evaluate (or approximate) 
the basic operations, namely the Racah transform and its associated 6j- 
symbol for arbitrary entries on the one hand, and the elementary braiding 
operators on the other. 

As a final remark we notice that the field theoretical approach gives us 
automatically the rate of growth of the absolute value of extended Jones 
polynomial with respect to the Chern-Simons coupling constant k = i + 2. 
The absolute value of the Reshetikhin-Turaev ^1] quantum invariants of 3- 
manifolds, [^(.M 3 )] (which are linear combinations of colored polynomials 
associated with surgery framed links C A4 3 ) are estimated to grow as 0(k d ), 
where the exponent d is bounded from above by some simple function (de- 
pending on the gauge group G) of the Heegaard genus of the manifold (cfr. 
E3|, Ch. 7). 
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Appendix 

A simple application of the procedure described in section 6 for the eval- 
uation of quantum link invariants is the explicit computation of the Jones 
polynomial for the plat presentation of the trefoil knot K tre { depicted in Fig. 
[*m The four strands are labelled by a same j\, together with its opposite j*, 
from left to right. 

Accordingly, the initial and final states to be associated with the quantum 
automaton are 



I[(((Ji,jT),(Ji,Ji)))];[0,0;-];0,0) 



(42) 



and 



|[(((jT,Ji),(Ji,Ji)))];[0,0;-];0,0) 



(43) 
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Figure 11: Plat presentation of the oriented trefoil knot. 



respectively, and they comply with the prescription of being odd vectors, see 
(|3HJ). From the picture we easily recognize that the operator to be employed 
is -Btref = (B 2 ) 3 . Since we are interested in the evaluation of the Jones 
polynomial we set from now on = |. Moreover, in order to apply 

the even braiding operator B 2 we have to perform preliminarly a duality 
transformation (|31|) on the odd vector 



|[y,i§Mo;-l;oo> = E 



1=0 



1 1 

2 2 

1 1 

2 2 



rl 1 1 1 
b' 2' 2' 2 



];[/;-]; 00), (44) 



which converts the initial state of the automaton into eigenvectors of the 
braid generator B 2 . The application of B tre f gives 



(S 2 ) 3 |[|,i,i,i];[0;-];00) = ^(A 



(+)\3 



1=0 



1 1 

2 2 

1 1 

2 2 



[2' 2> 2' 2]' ]' 



(45) 

where there appears the cube of the eigenvalue \\ + ^ defined in (jSHJ). Accord- 
ing to the expression of the extended Jones polynomial given in (|37|) and 
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taking into account (03J), we get 



V x (K tIei ; q) = [2] q ( [|, i, i, §]; [0; -]; 00 1 {B 2 f | [§, |, I ±]; [0; -]; 00> , 

(46) 

which, by using the orthogonality relations of the duality matrices (KAUL), 
amounts to 



V , (K tiei ; q ) = [2] q £ A[ +)3 = 1+ J + g3 (47) 

as required. 



Q 



37 



